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Abstract 

Coronal holes are the darkest and least active regions of the Sun, as observed both on the solar disk 
and above the solar limb. Coronal holes are associated with rapidly expanding open magnetic fields and 
the acceleration of the high-speed solar wind. This paper reviews measurements of the plasma properties 
in coronal holes and how these measurements are used to reveal details about the physical processes that 
heat the solar corona and accelerate the solar wind. It is still unknown to what extent the solar wind is fed 
by flux tubes that remain open (and are energized by footpoint-driven wave-like fluctuations), and to what 
extent much of the mass and energy is input intermittently from closed loops into the open-field regions. 
Evidence for both paradigms is summarized in this paper. Special emphasis is also given to spectro- 
scopic and coronagraphic measurements that allow the highly dynamic non-equilibrium evolution of the 
plasma to be followed as the asymptotic conditions in interplanetary space are established in the extended 
corona. For example, the importance of kinetic plasma physics and turbulence in coronal holes has been 
affirmed by surprising measurements from the UVCS instrument on SOHO that heavy ions are heated to 
hundreds of times the temperatures of protons and electrons. These observations point to specific kinds 
of coUisionless Alfven wave damping (i.e., ion cyclotron resonance), but complete theoretical models do 
not yet exist. Despite our incomplete knowledge of the complex multi-scale plasma physics, however, 
much progress has been made toward the goal of understanding the mechanisms ultimately responsible 
for producing the observed properties of coronal holes. 

1 Introduction 

Coronal holes are regions of low density plasma on the Sun that have magnetic fields opening freely into the 
heliosphere. Because of their low density, coronal holes tend to be the regions of the outer solar atmosphere 
that are most prone to behaving as a coUisionless plasma. Ionized atoms and electrons flow along the open 
magnetic fields to form the highest-speed components of the solar wind. 

The term "coronal hole" has come to denote several phenomena that may not always refer to the same 
regions. First, the darkest patches on the solar surface, as measured in ultraviolet (UV) and X-ray radiation, 
are called coronal holes. Second, the term also applies to the lowest-intensity regions measured above the 
solar limb, seen either during a total solar eclipse or with an occulting coronagraph. Third, there is a more 
theoretical usage that equates coronal holes with all open-field footpoints of time-steady solar wind flows. 
There are good reasons why these three ideas should be related to one another, but the overlap between 
them is not complete. To avoid possible confusion, this paper will mainly use the first two observational 
definitions, with the third one being only partly applicable. 

During times of low solar activity, when the Sun's magnetic field is dominated by a rotationally-aligned 
dipole component, there are large coronal holes that cover the north and south polar caps of the Sun. In 
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more active periods of the solar cycle, coronal holes can exist at all solar latitudes, but they may only persist 
for several solar rotations before evolving into different magnetic configurations. 

Despite not being as visually spectacular as active regions, solar flares, or coronal mass ejections 
(CMEs), coronal holes are of abiding interest for (at least) three main reasons. 

1 . The extended corona and solar wind connected with coronal holes tends to exist in an ambient time- 
steady state, at least in comparison with other regions. This makes coronal holes a natural starting 
point for theoretical modeling, since it often makes sense to begin with the simplest regions before 
attempting to understand more complex and variable structures. 

2. Coronal hole plasma has the lowest density, which makes it an optimal testbed for studies of colli- 
sionless kinetic processes that are the ultimate dissipation mechanisms in many theories of coronal 
heating. Other regions tend to have higher densities and more rapid Coulomb collisions, and thus the 
unique signatures of the kinetic processes (in their particle velocity distributions) are not as straight- 
forward to measure as in coronal holes. 

3. Coronal holes and their associated high-speed wind streams are also responsible for a fraction of 
major geomagnetic storms at 1 AU. Corotating interaction regions (CIRs) form when fast and slow 
wind streams collide with one another, and the subsequent interaction between these structures and 
the Earth's magnetosphere can give rise to long-lasting fluxes of energetic electrons. 

This paper reviews measurements of the plasma properties of coronal holes and how these measure- 
ments have been used to put constraints on theoretical models of coronal heating and solar wind accel- 
eration. There have be e n several earl ie r reviews that have focused m a inly on the topic of coronal holes, 
including 'Zirkei^ ('197?'), 'Suess' ('1979'), 'Harvev and Sheelev Jd (ll_979'), 'Parker' ('l99l'), 'Kohl and C ranmer 



(1999), Hudson (2002), Cranmer (2002a), Ofman (2005|), |dejSna and Arge (2005), Jones (2005?), and 
IWand (l2009h . Interested readers are urged to survey these other reviews in order to fill in any gaps in 
topical coverage in the present paper 

The remainder of this paper is organized as follows. Section |2] gives a brief history of the discovery 
and early years of research on coronal holes. Section [3] summarizes the observations and derived plasma 
properties of "on-disk" coronal holes (i.e., primarily using the definition of holes as dark patches on the 
solar surface at UV and X-ray wavelengths). Section previews the measurements of "off-limb" coronal 
holes and describes our current knowledge of how these structures are linked to various kinds of solar wind 
streams measured in situ. Section |5] discusses a broad range of possible theoretical explanations for how 
the plasma in coronal holes is heated and how the solar wind in these regions is accelerated. Section |6] 
concludes this paper with a few words about how the study of coronal holes helps to improve our wider 
understanding of heliophysics, astrophysics, and plasma physics. 



2 Historical Overview 



Although the term "coronal hole" was not first used until the middle of the 20th century, people have 
reported the existence o f visible features associa t ed with the Sun 's corona - seen during total eclipses - 
for centuries (see, e.g., IWang and Siscq 5ll980l: IVaquerol l2003h . A popular astronomy book from the 
first decade of the 20th century ( Serviss , 1909h contained clear descriptions of coronal streamers, eruptive 
prominences, and polar plumes in coronal holes. The following description of the latter, from an eclipse in 
1900, conveys that early speculation may sometimes be prescient: 

"The sheaves of light emanating from the poles look precisely like the 'lines of force' surrounding 
the poles of a magnet. It will be noticed in this photograph that the corona appears to consist of two 
portions: one comprising the polar rays just spoken of and the other consisting of the broader, longer, 
and less-defined masses of light extending out from the equatorial and middle-latitude zones. Yet even in 
this more diffuse part of the phenomenon one can detect the presence of submerged curves bearing more 
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Figure 1 : Left: Adaptation o f a sketch of the qui et so l ar cor ona made by PiddingtonI ( 1972h . based on prior 
drawings (IWaldmeien, 1 1 9551) and photographs (IGoldlll955h of the 30 June 1954 eclipse. Right: Contrast- 
adjusted eclipse image taken with the POISE instrument on 26 February 1998, in Westpunt, Cura9ao. The 
original image was made available courtesy of the High Altitude Observatory (HAO), University Corpora- 
tion for Atmospheric Research (UCAR), Boulder, Colorado. UCAR is sponsored by the National Science 
Foundation. 



or less resemblance to those about the poles. Just what part electricity or electro-magnetism plays in the 
mechanism of solar radiation it is impossible to say, but on the assumption that it is a very important part 
is based the hypothesi s that there exists a direct solar influence not only upon the magnetism, but upon the 



weather of the earth " dServ iss. 1909). 

The first quantitative observations of coronal holes were made by Waldmeier ( 19561 1957 ) at the Swiss 
Federal Observatory in Zurich. These features were identified as long-lived regions of negligible intensity 
in coronagraphic (off-limb) images of the 5303 A green emission line (see also Waldmeier, 1975, 1981). 
Waldmeier called the features that appeared more-or-less circular when projected onto the solar disk Locher 
(holes), and the more elongated features were called Kanal (channels) or Rinne (grooves). 

In off-limb eclipse and coronagraph images, the darkest coronal hole regions are surrounded by brighter 
and more complex streamers. These wispy structures appear to be connected to closed magnetic field lines 
at the solar surface, but they are often stretched upwards to an elongated cusp-like point, with thin "stalks" 
of radial rays at the top. For this reason their appearance was likened to a pointed German helmet (or a 
brush-topped Greek or Roman helmet), and the common phrase helmet streamers is often seen. The earliest 
studies of coronal morphology tended to concentrate more on streamers than cor onal h oles becau se the 
former are signi fi cantly easier to see than the latter (see, e.g.. Miller, 1908; Mitchell 1 1932 : Newkirk. ll967 



Pneumanlll968r) . iPiddington (Il972h outlined some early ideas about the global structure of the "quiet" (i.e.. 



solar minimum) corona. Figure [T] compares an adaptation of J. H. Piddington's sketch of the quiet corona 
to a more recent photograph from another eclipse around solar minimum. 

As the quality of the observations improved, coronal holes became objects of study in their own right. 
The largest coronal holes were observed to contain fine thread-like polar plumes that app ear to follow the 
superradially expand i ng open magnetic field lines above the solar limb (Saito, 1958; Stodd ard et al. . 119661: 



Newkirk and Harvey , 19681) . These e longated structures were found to correlate with bright chromospheric 



faculae on the surface (e.g.. iHarvevl 1 1 965 ) and wit h longer exten s ions f or the small jet-like spicules that 
continually rise and fall above the hmb jLippincott. 19571 : Beckersl 1968 ). 

Coronal holes were essentially re-di scovered in the l ate 1960s and early 1970s as discrete dark patches 
on the X-ray and ultraviolet sola r disk. Newkirk ( 19670 reviewed some of the earliest rocket-based mea- 
surements in the extreme UV, and lTousev et all ( Il968h discussed how the UV emission was "usually weaker 
over the poles" in images from a series of rocket flights between 1963 and 1967 (around solar minimum). 
These regions on the solar disk came to be called coronal holes in parallel with the earlier off-limb usage. 
Munro and Withbro3 (Il972) analyzed OSO-4 observations to conclude that both the density and electron 
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Figure 2: X-ray corona (0.25^.0 keV) observed by the Soft X-ray Telescope (SXT) on Yohkoh, on 6 
December 2000. Yohkoh is a mission of the Institute of Space and Astronautical Sciences in Japan, with 
participation from the U.S. and U.K. 



temperature were lower in these dark regions. In 1973 and 1974, solar instruments on the Apollo telescope 
mount (ATM) on Sk ylab confirmed m any earlier ideas about coronal holes with data significantly better in 
quantity and quality dHuber ef a/.L [l974; Kahler, 2000). 

In addition to the large north and south polar holes, there were also found to be smaller coronal holes 
that exist at lower latitudes (often at times other than solar minimu m). Sometimes the l a rgest coronal holes 
can exhibit thin "peninsulas" that jut out from the main regions. iHarvev and Recelv (2002) called these 
regio ns "polar lob es." Notable examples have been the so-called "Bo ot of Italy" seen by Skylab in 1 974 
(e.g.. lZirkeilll977l) and the "Elephant's Trunk" seen by SOHO in 1996 dOel Zanna and Bromage[ll999h . A 
third example, from December 2000, is shown in Figure ID 

Additional insights came from the fusion of spectroscopy and coronagraphic occultation. Inspired by 



rocket-borne UV observations of the extended corona during a solar eclipse in March 1970, IKohl et al. 



(Il978h developed a UV coronagraph spectrometer to measure the profile shape of the bright H I Lyman 
a emission line at 1216 A. This line is sensitive to several key properties of the velocity distribution of 
coronal protons, and thus these observations could be used to begin distinguishing proton temperatures 
from electron temperatures in the collisionless outer regions of coronal holes (see Section l43] ). The rocket- 
borne UV coronagraph spectrometer was launched three times (in 1979, 1980, and 1982), an d the results 



included the first dire ct evidence for proton h eating and supersonic outflow in coronal holes dKohl et al. 
1980l:IStrachan et ad lT993: Koh l ef a/.LE)06h . 



The fact that coronal holes coincide with regions of open magnetic field that expands out into inter- 
planetary space was realized during the fir s t decade of i n situ solar wind observations (e.g.. lWilcoxill968 : 



Altschuler et a/.LI1972l : lHundhausenlll972r) . iNocil d 19731) made a theoretical argument, on the basis of mea- 



sured wave fluxes and heat con duction, that coro nal holes should have the largest solar wind kinetic energy 
fluxes (i.e., the highest speeds). Pneumanl l 1973 ) argued that coronal holes need not have lower energy de- 
position than closed-field regions ( as is suggested by the lower intensities of coronal holes) if the solar wind 
carries away much of that energy. iKrieger et al\ dl973h utilized X-ray sounding rocket images to identify 
a large coronal hole as the solar source of a strong high-speed stream as measured by the Pioneer 6 and 
Vela spacecraft. Around this time it was also realized that coronal hole s and high-speed win d streams are 
also responsible for a fraction of the geomagnetic storms seen at 1 AU (Bell and Noci, 1973; Neupert and 
Pizzo" Tl974 : .Bell and Noel 197 6). (See also Tanskanen et al. C2005.) and Zhang et al. (.2007X ) Although 
there is still no complete understanding of which types of solar wind flow are connected with which types 
of coronal structures, the causal link between the largest coronal holes and high-speed streams remains firm 
(see also Section l^Tt . 



3 Properties of On-Disk Coronal Holes 

The traditional observational distinction between coronal holes and their surroundings (i.e., active regions 
and "quiet Sun") is that coronal holes have the lowest emission in the UV and X-ray. This definition must 
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be amended, however, to exclude filaments (which are often dark w hen projected again st the solar disk) 



that are c ool magnetic structures and not part of the corona (see, e.g.. Scholl and Habball 2008: Krista and 



lagr 

Gallagher. I2009h . Coronal hole magnetic fields are known to be more uniform and unipolar than in other 
regions (see below). The boundaries between coronal holes and surrounding re gions are sometimes sharp. 



sometimes diffuse, and sometimes filled with many small loops ( Hudson! 2002). 

Coronal holes are more or less indistinguishable from their surroundings in the photosphere and low 
chromosphere, and usually one cannot see any significant intensity contrast between hole and non-hole re- 
gions until the temperatu re exceeds about 10^ K. Spectra, however, can help make the distinction clearer 
Tephtskava et al. found that central self-reversals in the chromospheric Ca II H and K lines are 

noticeably different in coronal holes as opposed to surrounding quiet-Sun regions. A frequently used ob- 
servational diagnost i c of coronal hole boundar ies is the He I 10830 A near-infrared absorption line triplet 



(IHarvey et a/.LI1975tlHarvey and RecelyLl2002h . At these wavelengths, the absorption is weakest in coronal 



holes (i.e., the intensity is highest) and spectroheliogram images show the coronal holes quite clearly. It is 
somewhat counterintuitive that a spectral line of a neutral species (He") should be sensitive to the properties 
of the hot corona. However, the atomic level populations that determine the He I 10830 A source function 
are unusually responsive to photoionization from UV wavelengths shortward of 504 A. The overlying solar 
corona emits these wavelengths in abundance. In coronal hole regions, though, the corona generally has a 
lower density and temperature, and thu s there is less intense UV and X-ray emission to popul ate the upper 



levels of the He triplet state (see, e.g.. lZirinlll975l : lAndretta and Joneslll997l : ICenteno et al.luOOS) . This 



gives rise to reduced absorption. The He I 10830 A lines are also good probes of supersonic outflow veloc 



ities i n distant stellar winds (see, e.g.. lDupree et al. 
2008h . 



2005t iKwan et ail 120071 : ISanz-Forcada and Dupree , 



Another observational diagnostic of coronal holes is their elemental abundance signature (e.g.. Feld- 



man, 19981 : Feldman and Widingll2003 ). In the upper chromosphere, transition region, and low corona, 
holes exhibit abundances very close to those measured in the photosphere. This stands in contrast to other 
regions (quiet Sun and active regions), which show significant enhancements in the abundances of elements 
with low first ionization potential (PIP). This selective fractionation is believed to occur in the upper chro- 
mosphere, where low-FIP elements become ionized and high-FIP elements remain more neutral. These 
patterns extend into the heliosphere, where high-speed flows associated with coronal holes are often iden - 
tifiable from their near-photospheric FIP abundances (" von Steiger et al . 1995 : Zurbuchen et al . 2002). 
However, there is still no widespread agreement about the exact physical processes that give rise to this 
preferential ionization. 

The number, sizes, and heliographic locations of coronal holes vary as a function of the solar activity 
cycle. Large polar holes exist for about 7 years around solar minimum, and are not present for about 3 
or 4 years around solar maximum. However, in the declining phase of activity soon after the maximum, 
it is possible to see the gradual growth of the new-polarity polar coronal holes. This occurs as a number 
of sm aller high-latitude holes eventuallv collect together at the poles (.Timothv et al\. 19751: Harvev and 



tnsn 

Recely. l2002h . Taking this growth phase into account, there are only about 1 or 2 years at solar maximum 
without any distinct high-latitude coronal hole presence. Figure [3] illustrates the growth phase around the 
peak of solar cyc le 23 in 2001. The gr owth and full "reappearance " of polar coronal holes at this time 
was described bv lMiralles eTaZI (l2001bh and lMcComas et al l (l2002h . The post-maximum growth of new 
polar coronal holes lasts about twice as long as their disappearance in the rising phase of the next maximum 
(Waldmeier. 1981; Fish er and Sime , 1984). 

Many low-latitude coronal holes tend to be situated near the edges of magnetically complex active 
regions. Sometimes active regions emerge within the coronal holes themselves: these have been called "se a 
anemone" type regions from their spiky, flower-like appearance (Shibata ef al, 1994: Asai et a/.Ll2008h . 



Evolving magnetic interactions between active regions and coronal holes have been stu died both as a means 



of enhancing the mass flux of the associated solar wind on nearby open flux tubes (e.g., Habbal et al, 2008; 
Wang ef al, 2009) and as a p ossible explanation for the nearly rigid rotation of coronal holes (Wang et al , 
19961 : Antiochos et al . 2007). A burst of emerging magnetic flux in one of these active regions may give 



5 



Jan. 2001 Feb. 2001 Mar. 2001 Apr. 2001 




May. 2001 Jun. 2001 JuL 2001 Aug. 2001 



Figure 3: Polar view of the development of the north polar coronal hole from January to August 2001 (e.g., 
Carrington rotations 1972 to 1979), using reconstructed coronal hole boundaries from Kitt Peak He 1 10830 
A maps. The maximum of solar activity occurred between late 2000 and early 2001 . Data from the National 
Solar Observatory /Kitt Peak were produced cooperatively by NSF/NOAO, NASA/GSFC, and NOAA/SEL. 



rise to new systems of loop connections in the area bordering the coronal hole, and thus cause the coronal 
hole to decrease in size. This kind of rapid topological evolution of the magnetic field may be relevant in 
explaining extreme space weather events such as "t he day the solar wind disappeared" (i.e., the dramatic 



drop in the in situ density seen on 1 1 May 1999; see lJanardhan et a/.Ll2008h . 



Photospheric magnetograms show that coronal holes are more unipolar than other reg i ons; i.e., the y 
have a larger degree of magnetic flux imbalance between the two polarities ( Levine , 1982; Wand, 2009 ). 
For the large polar coronal holes, this appears to be the long-term outcome of the decay of active-region 
magnetic fields and their eventual diffusion up to the poles. The unipolar nature of coronal holes is likely 
to be related to their connection with open-field solar wind streams. As described in Section|2] one reason 
why coronal holes are dar k is that the solar wind car ries a way both mass and ene rgy, leaving a lower density 
and pressure. In addition. lAbramenko et al.\ (l2006ah and lHagenaar et al.\ (l2008h found that the local rate of 
emergence of small-scale magnetic flux (mostly in the form of "ephemeral" small-scale bipoles) is substan- 
tially lower in unipolar regions than in more mixed or balanced regions of positive and negative magnetic 
polarity. In most theories of coronal heating of closed loops, the total flux and the overall complexity of the 
field both drive the total heating rate. Thus, the lower emergence rate of new flux elements in coronal holes 
may be another factor determining why they have lower densities and pressures (i.e., less coronal heating; 
see, e.g.. 



Rosner et alX Il978f) and thus why they are dark. 



Over the past few decades, the magnetic connection between coronal holes and the fast solar wind has 
been traced down to ever smaller spatial scales. We now can say with some certainty that much of the 
plasma that eventually becomes the time-steady solar wind seems to originate in thin magnetic flux tubes 
(with observed sizes of order 50-200 km) observed mainly in the dark lanes between the ^1000 km size 
photospheric granulation cells. These strong-field (1-2 kG) flux tubes have been called G-band bright points 
and n etwork bright poin ts, and groups of t hem have b een sometimes ter med "solar filigree" (e.g. . Dunn and 
Zirker. 1 1 9731: ISpruia 1 1 984; Berger and TitleL l200ll; lTsuneta et ai, 2008). These flux fiibes are concentrated 
most densely in the supergranular network (i.e., the bright lanes between the larger ^30,000 km size su- 
pergranulation cells). Somewhere in the low chromosphere, the thin flux tubes expand laterally to the point 
where they merge with one another into a more-or-less homogeneous network field distribution of order 
~I00 G. This merged (mainly vertically oriented) field is associated mainly with the lanes and vertices be- 
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Figure 4: Summary of the largely unipolar magnetic field structure of polar coronal holes, with the fields of 
view successively widening from flux tubes in intergranular lanes (a), to a "funnel" rooted in a supergranular 
network lan e (b). and finallv to the extended corona (c). Adapted from Figure 1 of Cranmer and van 
Ballegooijen (l2005h . 



tween supergranular cells. Because this field does not fill the entire coronal volume, it is still susceptible to 
an additional stage of lateral expansion and weakening. Thus, at a larger height in the chromosphere, these 
network flux bundle s are believed to undergo further broadening into so-called "funnels' ' (iGabrieill976l: 
Dowdy et al. , However, it is still unclear to what e xtent the closed fi elds in the supergranular cell 



centers affect the canopy-like regions between funnels (e.g.. lSchriiver and T itle. 2003). Figure |4] illustrates 
the successive merging of unipolar field in coronal holes. 

Observations of blueward Doppler shifts in supergranular network lanes and vertices, especially in 
large coronal holes, may be evidence for either the "launching" of the so lar w ind itself or fo r upward- 
going waves that are linked to wind acce leration processes (e.g., iHassler et all 1 19991 : IPeter and Judg3, 



19991; lAiouaz et a/.L 120051 : iTu et a/.Ll2005h . These measurements are consistent with several models of the 



dynam ic, multi-species solar wind in superradially expanding funnels (Bvhr ing et al.i 2008; Marsch et al. 



2008r) . However, this interpretation of the data is still not definitive, because there are other observational 



diagnostics that have sh own more of a blueshif t in the sup ergranular cell-centers between funnels (e.g.. He 
I 10830 A: iDupree et a l. 1996; Malanushenk o and Jone s. 2004). There may be subtle radiative transfe r 



effects that preferentially brighten regions of upflow or downflow (see, e.g.. lChae etal. , 1997 ; Avrett , 1999h 



and these may need to be taken into account in order to understand the meaning of the measured Doppler 
shifts in these regions. 

Lastly, it is important to mention the phenomenon of transient coronal holes (sometimes known as 
"coronal dimmings"). These are rapid reductions in the UV and X-ray intensity in discrete patches that 
appe ar to b e spati ally and temporally correlated with the hftoff of CME plasma (e .g.. Rust, 1983; Thompson 
etal. ]2000l:lYange/adl2009h . Spectroscopically, t hese regions are see n to exhibit similar characteristics as 



normal coronal holes, including Doppler blueshifts (iHarra et alX . 120071) and large ampUtudes of nonthermal 



wave broadening (Mcintosh, 2009). UV coronal dimmings are beginning to be used as diagnos t ics for the 



amount of plasma released - i.e., the total mass - in the CME event (e.g.. lAschwanden et al. . l2009h . It 



should be noted that transient coronal holes represent just one kind of observed dimming that is associated 
with ti me-dependent flare/CME activity; there are other types of dimmings that do not resemble coronal 
holes (iHudsonll2002h . 
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4 Properties of Off-Limb Coronal Holes 



Coronal holes observed above the solar limb usually trace out the same regions that are identified as dark 
coronal-hole "patches" directly on the solar disk. Thus, the lower plasma density that more or less defines 
the off-limb coronal hole is directly related to the lower density measured on the disk. Section |4T|briefly 
discusses how these regions are believed to be connected magnetically with the broader heliosphere. The 
observations of off-limb coronal holes made with visible-light imaging and polarimetry (Section [4.2b and 
ultraviolet spectroscopy (Section l43T l are also summarized below. 



4.1 Magnetic connectivity with the solar wind 



Although the magnetic field in the solar corona is generally too weak to be measured directly, the overall 
morphology of the field lines can be extrapolated from magnetograms taken at the level of the photosphere. 
One popular technique is the "potential field source surface" (PFSS) method, which assumes the corona 
is current-free out to a spherical surface (set typically at a radius between 2.5 and 3.5 solar radii, or ^0), 
above wh ich the field is radial (e.g., Sch attenef g/., 196 9; Altschuler and Newkirk, 1969; Hoeksema and 
Scherrer, Th e PFSS method has been shown to create a relatively good mapp i ng be t ween the Su n 

and the heliosphere ("Arge and Pizzo','2000'; 'Luhmann et aZ.V2002l;IWang and Sheeley Jr,'2006';'W an^ l2009h. 
althou gh the results can b e problematic for regions dominat ed bv stream-stream interac tions (Poduval and 
Zhao^laOOj. 

By far, the strongest causal link between a specific type of coronal structure (measured via remote 
sensing) and a particular type of quasi-stead y solar wind flow (measured in situ ) is the connection between 
large coronal holes and high-speed streams dWilco xL ll968l:lKrieger et a/.LlT973h . The general interpretation 
of this correlation - together with the results of magnetic extrapolation models such as PFSS - is that coronal 
holes represent a bundle of open flux tubes that flare out horizontally as distance increases. In other words, 
the coronal hole flux tubes expand superradially. Although there are some observations that appear to 



suppo rt other interpretations dWoo et al. , ll999l ; lHabbal et adl200l l ;IWooll2005l;IWoo a nd Druc kmullerovl 



2008r) . the preponderance of evidence seems clearly to s upport the idea that fast solar wind streams emerge 
mainly froni superradially expanding coronal holes (e . g.. | Guhathakurta et al. ■ Il999al: l ICranmer et al.\ 



Jonesll2005l;[Wang and Sheeley Jill2006l; IWang ef a/.L l2007h . 

In contrast to the rather definitive correlation between large coronal holes and the fa st solar wind, th e 
coronal sources of the more chaotic slow-speed solar wind are not as well understood (see lSchwenn , l2006h . 
Two regions that are frequently cited as sources of slow wind are: (1) boundaries between coron al holes 



and s treamers, and (2) narrow plasma sheets that extend out from the tops of streamer cusps (Wang et al 
l200(]|;IStrachan et al 1 12OO2I; ISusino et ali ilOO^ . These regions tend to dominate around solar minimum. 
Note that the former type of boundary region tends to contain flux tubes that may be classified as coronal 
holes when using the theoretical definition (i.e., footpoints of field lines that are open; see Section[Tli but 
would not be defined as such when using the observational definitions (i.e., low emission or low density). 

During more active phases of the solar c ycle, th ere is evidence that slow solar wind streams also emanate 
from s mall coronal holes (e.g.. [Nolte et al. . 19761; N eugeb auer et al, 1998; Zhang et al, 2003) and active 



regions (IHick et a/.LIl995tlLiewer et a Sakao et 



2007h . During the rising phase of solar activity, 
there seems to be a relatively abrupt (< 6 month) change in the locations of slow wind footpoints: from 
the high-latitude hole/s treamer boundaries to the low-latitude active region and small coronal hole regions 



dLuhmann et a/.Ll2002h. The ability of many of these kinds of regions to produce slow wind was modeled 



by lCranmer et all (l2007h and lWang et al.\ (120091) ; see also Section|5] 



4.2 Visible light observations 

Measurements of the plasma properties in the extended corona (i.e., r w 1 .5 to 10 Rq, where the main solar 
wind acceleration occurs) require the bright solar disk to be occulted. The coronal emission is many orders 
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Figure 5: Comparison of empirically determined densities in the upper solar atmosphere. I Avrett and Loeser 
( l2008r) values of electron number density (solid black curve) and total hydrogen number density (dot-dashed 
black curve) are compared with various visible-light pB electron number densities for coronal holes (blue 
curves) and streamers (red curves); see text for details. 



of magnitude less bright than the emission from the solar disk, so even the vast majority of "stray light" that 
diffracts around the occulting edge must be eliminated. Visible-light coronagraphs that combine stray-light 
rejection with linear polarimetry have the ability to measure the Thomson-scattered polarization brightness 
(pB) in the corona. The use of pB rather than the total coronal brightness eliminates the contribution from 
the dust-scattered F-corona, which is believed to be unpolarized up to distances of about 5 Rq. Because 
the coronal plasma is optically thin to the Thomson-scattered photons, pB is proportional to the line-of- 
sight integral of the electron density n,,, multiplied by a known scattering function. Methods for inverting 
this integral to derive as a function of position in various coronal structures have been developed and 
improved over the years (e.g., van de Hulst','l950"; ' Altschuler and PerryLll972l : lMunro and JacksonLllQ??! : 
[Ouhathakurta and Holzer, 1994; Frazin et ai, 200^ For coronal holes, the LASCO (Large Angle and 
Spectrometric Coronagra ph) instrument on SOHO ha s also been used to probe the superradial expansion of 
open ma gnetic flux tubes ( DeForest etal . 1 997 , 200 1 ) and the evolution of transient polar jets ("Wang et al. , 
[l998; Wood ef adll999h . The White Light Coronagraphs on Spartan 201 (Fisher and Ouhathakurta, 1995) 
and on the UVCS ( Ultraviolet Coronagraph Spectrometer) instrument aboard SOHO (e.g.. lKohl et a/.LIl995 ; 
iRomoli ef ad 12002 ) have provided electron densities between L5 and 5 Rq in coronal holes. 

Figure |5] illustrates a selection of visible-light measurements of the electron density in coronal holes 
and compares them to similar measureme nts of streamers and to a semi-empirical model of the chromo- 
sphere, transition region, and low corona ('Avrett and Loesei*, '2008') . The blue cor onal hole curves were 
adap ted from the r esults of Fisher and Ouhathakurta ( 1995) (dotted), Cran mer et al\ ri 999b) (solid), Doyle 
et al. (Il999h (dashed), and lOuhathakurta e t al. ( 199 9bl) (dot-dashed ) . The red curves for equa torial helme t 
streamers were adapted from the results of ISittler and Ouhathakurta ( 1999h (solid) and lOibson et al. ( 1999h 
(dashed). 

Note that streamers are denser than coronal holes by about a factor of 10, but the hole measurements 
themselves can often exhibit variations in the electron density by factors of the order of 2-3. Much of 
this spread is likely to be the result of different lines of sight passing through regions that contain varying 
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numbers of polar plumes (see, e.g. JCranmer et al. . ll999hL l2008h . Some of the cited pB observations were 
optimized to avoid bright concentrations of plumes, and others have been purposefully averaged over the 
full range of coronal hole substructure. It is also possible that absolute calibration uncertainties may still 
exist between the different instruments used to determine pB and tie, and this could compound the reported 
range of variation in coronal hole electron densities. 

For a steady-state solar wind outflow, the conservation of mass demands that the product of density, flow 
speed, and cross-sectional area of the flux tube remain constant. Thus, if the magnetic geometry and the 
elect ron density are known, mass conservation allows the solar wind outflow speed to be computed. Kohl 
et al. (2006) used the representative values of shown in Figure |5] together with a range of estimates for 
the su perradial flux-tube expansion of coronal holes to determine outflow speeds in coronal holes. Figure 
41a of Kohl et al. ( 2006h illustrates the result of this process, which shows a large range of values reflecting 
the uncertainties in both and the flux-tube area factor. Despite these uncertainties, though, the electron 
densities that became available in the 1990s demonstrated that the fast solar wind accelerates rapidly in 
coronal holes - probably reaching half of its asymptotic terminal speed (Mqo ~ 700-800 km s"^) by heights 
no larger than 2-4 Rq. 

The increased sensitivity of the LASCO instrument over earlier coronagraphs revealed a nearly contin- 
ual release of low-contrast density inhomogeneities. or "blobs." from the cusps of helmet streamers ( Shee- 



ley Jr et a/.. ll997l : lTappin et a/.L Il999VWang e t ad boOOl: IChen et a/.Ll2009h . These features were seen to 
accelerate up to speeds of order 300-400 km s"' by the time they reached the outer edge of the LASCO field 
of view (r « SO/?©); see also Figure [8] below. The blobs are typically only about 10% to 15% brighter or 
dimmer than the surrounding streamer material. Because of this low contrast, these features do not appear 
to comprise a large fraction of the mass flux of the slow solar wind. However, it is still unclear whether 
blobs are passive "tracers" that flow with the solar wind speed, or whether they are wavelike fluctuations 
that propagate relative to the background solar wind reference frame. This diagnostic tool has been much 
more difficult to apply in coronal holes than it has in the bright streamers, so no firm measurements of the 
fast wind acceleration yet exist from this technique. 

Visible light measurements have also revealed evidence for compressive magnetohydrodynamic (MHD) 
waves that propagate along open field lines in coronal holes. Intensity oscillations measured by the UVCS 
and EIT instruments on SOHO were found to have periodicities between ab out 10 and 30 minutes and are 
consistent wi th being upwardly propagating slow-mode magnetosonic waves (iDeForest and Gurmanlll998 : 
Ofman et a/.Lil999. 2000). The relative amplitu de of the density fluctuations (5n/no) for these waves was 



found to range between about 0.03 and 0. 15 (see Cranmeil 2004a). This is consiste nt with measurem ents of 
the density fluctuation amplitudes made at larger distances via radio scintillations (Spanglen. l2002h and in 
situ instruments (Tu and Marsch. .L9 94). There have also been claims that low-frequency oscillations h ave 
been measured in H I Lya emission dMorgan et all l2004l : iBemporad et al.i l2008l : iTelloni et al. , l2009h . In 
these cases, however, it is extremely important to take into account all of the relevant instrumental effects. 
These measurements still appear to be provisional. 

As seen in Section |2] above, coronal holes have long been observed as the sites of thin, ray-like polar 
plumes. The earliest measurements of polar plume properties were made in broad-band visible light, and 
these dense strands are often seen to stand out distinctly from the ambient interplume corona. It is not clear, 
though, to what extent off-limb observations (which integrate over long optically thin lines of sight) ever 
capture only the "pure" plume or interplume plasmas. Space-based observations from, e.g.. Spartan 201 
and SOHO i mproved our abilitv to measure the physical properties in and between plumes (e.g.. Fisher and 
Guhafliakurta. ll995l:lDeForest ef a/.LIl997l:ICranmer ef a/.LIl999bl : lDeForest et a/.Ll200lh . Although the 
brightest plumes are still discernible at the uppermost heights observed by LASCO (i.e., 30^0 Rq), the 
plume/interplume density contrast becomes too low to measure clearly in interplanetary space (r > 6QRq). 
How ever, indirect and possibly plume-related signatures in the in situ data have been reported by Thieme 



et al. (Il990h .l Reisenfeld et al .l il999h . and lYamauchi et al\ (l2002h . The disappearance of plumes probably is 
the result of some combi nation of transver s e pres sure balance (i.e., dense plumes expanding to fill more of 
the available volume; see lDel Zanna et adll998h and MHD instabilities that can mix the two components 
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dParhi et a/.[ll999l:lAndries and Goossensll200lh . 



4.3 Ultraviolet spectroscopy 

Ultraviolet spectroscopy of the c orona is a powerful tool for obtaining detailed em pirical descriptions of 



solar plasma conditions (see, e.g.. lKohl and Withbroelll982tlWithbroe et a/.LI1982h . Coronal holes, being 



the lowest density regions of the outer solar atmosphere, exhibit a complex array of plasma parameters due 
to their nearly colhsionless nature. As a result, every ion species tends to have its own unique temperature, 
its own type of departure from a Maxwellian velocity distribution, and its own outflow speed. After a brief 
summary of near-limb measurements made with the SUMER instrument on SOHO, this section mainly 
describes results at larger heights (in the more clearly colhsionless extended corona) from UVCS. 

The un-occulted SUMER (Solar Ultraviolet Measurements of Emitted Radiat ion) spect rometer has ob- 
served off^imb regions up to heights of approximately 1.3 Rq in coronal holes ( Wilhelm et al. 1995.2OO0l 



20041 120071) . In polar regions at solar minimum, ion temperatures exceed electron temperatures even at 
r ~ \ .05Rq, where densities were presumed to be so high as t o ensure rapi d coUisional c oupling and thus 



equal temperatures for all species (ISeelv et alX. 119971 : iTu et al.i 11998: Landi and Cranmer , 12009) . Spectro 
scopic evidence is also mounting for the presence of transverse Alfven waves pro pagating into the corona 
dBanerjee et aIlll998[l2009l: lDolla and Sol omon'. '2008': 'Landi and Cranmer. '2009|)^ 



E lectron temperatures derived from line ratios (Habbal et ai. 1993 : David et alx Il998c iDoschek et al. 



200ll:IWilhelriiLl2006l : lLandiL l2008h exhibit relatively low values in the off-limb coronal holes {Tg ~ 800, 000 



K) that are not in agreement with higher temperat ures derived from "frozen-in" in si tu charge states (Ko 
et q/.. .1997) . It is difficult to reconcile these observations with one another in the absence of either non- 
MaxweU ian electron velocity di stributions or strong differential flows between different ion species near 
the Sun (lEsser and EdgailEoOol) . However, some m odels consiste nt with both the SUMER temperatures 
and the in situ charge states are being produced (e.g.. lLaming and Le pri. 2007). 

Figure|6]displays a range of temperatures measured in coronal holes and high-speed wind streams, and 
it shows how the SUMER electron temperatures {T^) are noticeably lower than the heavy ion temperatures 
(e.g., shown for the O"^^ ions that correspond to the bright O VI 1032, 1037 A spectral lines) even in the 
low corona. The degree of agreement between th e spectroscopic measurem ents and one-fluid models of 
the low corona (i.e., the semi-empirical model of ' Avrett and Loesei (2008) and the theoretical model of 



Cranmer et al\ ( 120071) ) depends on the height of the sharp transition region between chromospheric (10"* K) 
and coronal (10^ K) temperatures. 

The UVCS instrument on SOHO is a combination of an ultraviolet spectrometer and a linearly occulted 
coronagraph that observes a 2.5 Rq long swath of the extended corona, or iented tangentially to the radial 



direction at heliocentric distances ranging between about 1.5 and 10 Rq (IKohl et a/.LI1995L 119971 11998 



20061) . In coronal holes, UVCS measurements have allowed many key details about the velocity distribu- 
tions of H'\ O^^, and Mg^' to be derived. For the resonantly scattered emission lines seen at large heights 
with UVCS, the most straightforward plasma diagnostic is to use the Doppler-broadened line width as a 
sensitive probe of the overall variance of random particle motions along the line of sight. In other words, 
measuring the line width provides a constraint on the so-called "kinetic temperature" (i.e., a combination 
of microscopic stochastic motions and macroscopic [but unresolved] motions due to waves or turbulence) 
along the direction perpendicular to the (nearly radial) magnetic field lines. 

In the ionized solar corona, a given hydrogen nucleus spends most of its time as a free proton, and only 
a small fraction of time as a bound H*^ atom. Thus, the measured plasma p roperties of neutra l hydrogen are 
considered to be valid proxies of the proton properties below about 3 Rq (lAllen et al. [ 12000^. Spartan 201 
and UVCS observations of the H I Lya emission line in coronal holes indicated rather large proton kinetic 
temperatures in the direction perpendicular to the magnetic field {Tp±^ ~ 3 MK) and also the possibility of a 



1997: Cranmer et al. 



mild temperature anisotropy (wit h T„±^ > Tm) above heights of 2-3 Rq dKohl et al 
[l999t; Antonucci et al... 2004: Kohl e7a/.l 2006). 

UVCS observations indicated that the O"^^ ions are much more strongly heated than protons in coronal 



11 




Figure 6: Radial dependence of empirical and model temperatures in polar coronal holes and fast wind 
stream s. Mean plasma temperatures from a semi-empirical model (dashed blac k curve; Avrett and Loesen , 
2008 ) and from a turbulence-driven coronal he ating mod el (solid black curve; Cranrner et al\. 2007 ). Te 



from off-limb SUMER measurements made by Wilhelm (2006) (dark blue bars) and Landil ( 2008h (light 
blue bars), 7), from UV CS measurements as s emble d by lCranmer (,2004b.) (se e text), and perpendicu lar O"*"^ 
ion temperatures from lLandi and Cranmer (l2009h (open green circles) and ICranmer et al. \ (l2008h (filled 



green circles). In situ proton and electron temperatures in the fast wind (r > 60Rq) are from lCranmer et al. 

Eqq3). 



holes, with perpendicular temperatures in excess of 200 MK (see Figure |7|. This exceeds the temperature 
at the central core of the Sun by an order of magnitude! The UVCS mea surements also provided signa - 
tures of temperature anisotropics possibly greater than T±i/T\\i « 10 (e.g.. lCranmer et al. I ll999bLE)08h . 
The measured kinetic temperatures of O"*"^ and Mg"^^ are significantly greater than mass-proportional when 



compared with protons, with Tj/Tp > mi/nip (see also lKohl et a/.LIl999ll2006r) . The surprisingly "extreme" 



properties of heavy ions in coronal holes have led theorists to develop a number of new ideas regarding the 
heating and acceleration of the solar wind; these are discussed further in Section l574l 

Figure |6] shows UVCS perpendicular temperatures for protons and O"^^ ions in coronal holes. The O^^ 
data points were taken from the recent re-analvsis of solar m inimum data from 1996-1997 bv Cranmer 
et al. (l2008h . The proton temperature data were assembled by Cranmer (2004b|) from a number of indi- 
vidual measurements o f the H I Lyg profile a t solar minimum. T he sources of these m easurements are: 
Cranm er et al. ' ('1999b') (squares), lEsser et al. 1(11999) (diamonds), IZangrilli et a /1 (ll999h (asterisks), and 
Antpnucci et al. (2000) (triangles). The kinetic proton temperatures are of order 2-3.5 MK, but in Figure|6] 
we attempted to remove the co ntribution of nonthermal wave broadening. The semi-empirical model of 



Cranmer and van BallegooiienI (12005 ) was used to specify the amplitude of transverse Alfven waves as a 
function of height, and their contribution to the line widths was subtracted. The remaining "microscopic" 
T±p does not show as clear a signal of "preferential" proton heating as would be apparent from the larger 
kinetic temperature. Although one can still marginally see that Tp > Tg, the existing measurements of Tp 
and Te do not fully overlap with one another in radius. Improved measurements are needed in order to better 
constrain the proton and electron heating rates in the corona. 

The UVCS emission line data contain information about the Doppler motions of atoms and ions along 
the magnetic field (i.e., transverse to the line of sight). The so-called "Doppler dimming" diagnostic tech- 



12 






1 / . 






1 J . 






lff25 1030 1035 1040 
wavelength (A) 



1045 1025 



1030 1035 1040 
wavelengtb (A) 



1045 



On-disk profiles: T = 1-3 MK 



Off-limb profiles: T ^ 200 MK 



Figure 7: Combined image of the solar corona from 17 August 1996, showing the solar disk in Fe XII 195 
A intensity from BIT (yellow inner image) and the extended corona in O VI 1032 A intensity from UVCS 
(red o uter image). Axisymmetric field lines are from the solar-min imum model o f Banaszkiewicz et al. 
1 1998V and O VI emission line profiles (bottom) are from SUMER dWarren et a/.LIl99l left) and UVCS 
(iKohleraniT997[ right). 



nique provides constraint s on both the bulk outflow speed along the field and the parall el kinetic tempera- 



tie parall el 

ture (for more details, see iKohl and Withbroelll982l : lNoci et ami987tlKohl et oZlbooeh . In coronal holes, 
Doppler-dimmed line intensities from UVCS are consistent with the outflow velocity for O^^ be i ng larger 
than the outflow velocitv for protons bv as much as a factor of two at large heig hts (Ko hl et a/.L 1998: Li 
et a/.. ll998l:ICranmer et adll999bh . Figure [8] illustrates the outflow speeds measur ed by UVCS in corona l 
holes, and compares with the theoretical model of the fast solar wind presented bv lCranmer et al\ (120070 . 
Also shown for comparison are observational and theoretical data for the slow solar wind associated with 
equatorial helmet streamers at solar minimum. 

In contrast to many prior analyses of UVCS data, which concluded th at there must be both intense 
preferential heating of the O^^ i ons and a strong fiel d-aligned anisotropy, iRaouafi and Solankil (120041) . 



Raouafi and Solanki (l2006h . and iRaouafi et al\ ( 120071) reported that there may not be a compelling need 



for O"""^ anisotropy depending on the as sumptions made about the other plasma properties of the coronal 
hole (e.g., electron density). However, ICranmer et al\ (l2008h performed a detailed re-analysis of these 
observations and concluded that there remains strong evidence in favor of both preferential O"^^ heating 
and acceleration and significant O^^ ion anisotropy (in the sense Tu > T\\i) above r « 2.1 Rq in coronal 
holes. In determining these properties, it was found to be important to search the full range of possible 
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Figure 8: R adial dependence of solar wind outflo w speeds. UVCS Do ppler dimming determinations for 
protons (red: lKohlef g/. . l2006h and O ions (green: ICranmer et a/.Ll2008h are shown for polar coronal holes, 
and are compared with theoretical models of the polar and equatorial solar wind at solar minimum (black 
curves; ^ranmer et ai, 20 07|) and th e speeds of "blobs" measured by LASCO above equatorial streamers 
(open circles: iSheeley Jr et al.i\l99m . 



ion temperatures and flow speeds, and not to make arbitrary assumptions about any given subset of the 
parameters. 

The UVCS results discussed above are similar in character to in situ measurements made in the fast so- 
lar wind, but they imply more extreme departures from thermodynamic equilibrium in the extended corona. 
For example, proton velocity distributions measured in the fast solar wind between 0.3 and 1 AU have 
anisotropic cores with Tp± > r^,[], and their magnetic moment s increase with increasing dis t ance; this im- 
plies net input of perpendicular energy on kinetic scales (e.g.,'Marsch et a/.Vl982'; Marsch,'2006V Many 
heavy io ns appear to flow faster than the protons b y abo ut the local Alf y en spe ed (Hefti etai, 1998; Reisen- 
feld et al., 1200 ih . and in the fastest s olar wind streams they are also heated preferentially in the same sense 
as in the corona I CoUier eTal. , [1996'). 

In the years since the solar minimum of 1996-1997, UVCS observed a large number of other coronal 
holes that appeared throughout the maximum of solar cycle 23 and the new-millennium solar minimum of 
2007-2009. UVCS tends to observe only the largest coronal holes, since when the smallest holes rotate to 
the solar limb their UV line profiles tend to be contaminated by emission from streamers in the foreground 
and background. This selection effect naturally screens out small coronal holes that have been correlated 
with slow solar wind streams at 1 AU (' Nolte et 'al. , 1 19761 : iNeugebauer et ail 119981) . In the cases where 
UVCS and in situ measurements were made for the same regions associated with large coronal holes, hig h 



speeds in excess of 600 km s were inevitably seen in interplanetary space (iMiralles et all l2004l 1200^ 



However, the O"^^ outflow speeds measured in the extended corona via Doppler dimming have showed a 
substantial range of variation. For example, Miralles et al. (2001a) found that the outflow speeds at 2-3 
Rq in an equatorial coronal hole were approximately thr ee times lower than those measured in the polar 
coronal holes from 1996-1997 at similar heights (see also lPoletto et a/.M2002h . This implies that the range 
of coronal heights over which the acceleration of the solar wind occurs can vary greatly, even when the 
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wind at 1 AU ends up similarly fast. 

UVCS also has measured the plasma properties in bright polar plumes. The densest concentrations of 
polar plumes along the line of sight are seen to ex hibit narrower line widths - i.e., lower kinetic tempera - 
tures - than the lower-density interplume regions ( Kohl et al. , 1997 ; Noci et al. , 19971: Kohl et a/.L 12006 ). 
Similarly, plumes are seen to have lower outflow speeds than the interplume regions Joiordano et all 2000L 
iTeri aca et al.\ 2003; Raouafi et al. 2007), althouah at low heights the data a re not as clear-cut fe.g.. Gabriel 
et al.,2£>03). UVCS also has put c onstraints on the nlume/i nterplume density contrast and the filling factor 
of polar plumes in coronal holes. ICranmer et all (Il999g) used a large number of UVCS synoptic mea- 
surements to determine statistically that, at heights around r « 2Rq, the plume/interplume density ratio is 
approximately 2, and polar coronal holes are comprised of about 25% plume and 75% interplume plasma 
(corresponding to ~40 individual plumes distributed throughout the coronal hole). Earlier measurements 
made closer to the limb showed a higher density contrast and a smaller filling factor, so the UVCS data are 
generally consistent with lateral expansion of polar plumes with increasing distance. 

It is still relatively unknown how much of the mass, momentum, and energy flux of the fast solar wind 
comes from polar plumes. Despit e that uncerta i nty, t hough, there have been several reasonably successful 
models of polar plume formation. Wangl l 1994 . 1998i) presented models of polar plumes as the extensions 
of flux tubes with concentrated bursts of added cor onal heating at the base - presumably via nanoflare-like 
reconnection events in X-ray bright points (see also^ DeForest et al. , 2001 ). In these models, the extra basal 
heat input is balanced by conductive losses to produce the larger plume density. The heating rate in the 
extended corona is not affected by the basal burst, but the larger density in the flux tube gives rise to less 
heating per particle at all heights, which leads to lower temperatures in the extended corona and a smaller 
gas pressure force for solar wind acceleration. This model is consistent with the smaller temperatures and 

outflow speeds measured in plumes with UV spectroscopy. 

_UVCS made the first spectroscopic measurements of polar jets in coronal holes ("Pobrzvcka et al. , 2000l 



2006) . The events observed by EIT, LASCO, and UVCS during the 1996-1997 solar minimum tended to be 



"cool jet s" with higher dens ities, lower temperatures, and faster outflows than the surrounding coronal holes 
(see also Wang et flZ.,'l998'). More recently, Hinode has obse rved a ne w population of "hot" X-ra y jets in 
coronal holes (Culhane et al. 2007; Cirtain et al. 2007 ; Shimojo et a/.L [2007 ; Filippov ef a/.Ll2009h . UVCS 
found that some of these events per sist up to heights of at least 1 .7 Rq and that the jet protons remain hotter 
than the surrounding coronal hole (iMiralles et a/.ll2007n . Thus, there appear to exist two distinct kinds of 
polar jets (cool and hot), with differences possibly related to the relative degrees of heating and adiabatic 
expansion of the jet parcels. Jets and plumes have roughly similar angular sizes and i ntensity contrasts nea r 
the solar limb, and there is growing evidence that they share a common origin (e.g.. lRaouafief a/.L 12008*). 
It is possible that the only substantial difference between the two phenomena is the duration of the bursts 
of basal heating; i.e., jets seem to be the result of short-lived bursts of heating, whereas plumes may be the 
product of base-heating events that last longer than several hours. 



5 Coronal Heating and Solar Wind Acceleration 

Despite more than a half-century of study, the basic physical processes that are responsible for heating the 
million-degree corona and accelerating the supersonic solar wind are not known. This section broadens 
the topic of this paper a bit beyond just coronal holes, since an understanding of solar wind acceleration 
naturally encompasses not only the question of why fast solar wind streams are fast, but also why (various 
kinds of) slow solar wind streams are slow. Section ISTI summarizes some of the major issues regarding 
coronal energy deposition. The next two subsections describe two alternate views of solar wind acceleration 
via waves and turbulence in open flux tubes (Section |521) and reconnection between open and closed flux 
tubes (Section l5.3l l. Lastly, Section \5A\ reviews how collisionless kinetic effects in coronal holes (i.e., 
preferential ion heating and temperature anisotropics) can be used to more conclusively identify the detailed 
physical processes that produce the solar wind. 
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5.1 Sources of energy 

Different physical mechanisms for heating the corona probably govern active regions, clos ed loops in the 
quiet corona, and the open field lines that give rise to the solar wind (see other reviews by Marschl 1999t : 



HoUweg and Isenberd, 2002 ; Longcope , 2004 ; Gudiksen , 2005 ; Aschwan den. 2006; Klimchuk, 2006 )■ The 



ultim ate source of the energy is the solar convection zone (e.g.. lAbramenko et a/., ,2006b; Mcintosh et a. 



2007n . A key aspect of solving the "coronal heating problem" is thus to determine how a small fraction of 
that mechanical energy is transformed into magnetic free energy and thermal energy above the photosphere. 
It seems increasingly clear that loops in the low corona are heated by small-scale, intermittent magnetic 
reconnection that is driven by the continual stressing of their magnetic footpoints. However, the extent to 
which this kind of impulsive energy addition influences the acceleration of the solar wind is not yet known. 

Intertwined with the coronal heating problem is the heliophysical goal of being able to make accurate 
predictions of how both fast and slow solar wind streams are acce lerated. Empirical correlation tech niques 
have become more sophisticated and predictivelv powerful (e.g., Wang and Sheelev Jr, 1990, 2006; Arge 
and Pizzo. i2000:.Leamon and Mclntoslu, 2007, ; , Cohen et a/.„,2007;,Vrsnak et al.^.lQQl) but they are limited 
because they do not identify or utihze the physical processes actually responsible for solar wind accelera- 
tion. There seem to be two broad classes of physics-based models that attempt to self-consistently answer 
the question: "How are fast and slow wind streams heated and accelerated? " 

1 . In wave/turbulence-driven (WTD) models, it is generally assumed that the convection-drivenjostling 
of magnetic flux tubes in the photosphere drives wave-like fluctuations that propagate up into the ex- 
tended corona. These waves (usually Alfven waves) are often proposed to partially reflect back down 
toward the Sun, develop into strong MHD turbulence, and dissipate over a range of heights. These 
models also tend to explain the differences between fast and slow solar wind not by any major dif- 
ferences in the lower boundary conditions, but instead as an outcome of different rates of lateral 

I 1 1 ' 

flux-tube expansion over several solar radii as the wind accelerates (see, e.g., Hollweg, 1986; Wang 
and Sheeley Jr, J991 ; Mat t haeus et al., 1999; Cranmer, 2005; Su zuki, 2 006; .Suzuki and InutsukaL 
l2006l : E;ranmer ef aZ.[[2007l:IVerdini and Velli 120071: 1 Verdini ef alll2009h . 



2. In reconnection/loop-opening (RLO) models, the flux tubes feeding the solar wind are assumed to 
be influenced by impulsive bursts of mass, momentum, and energy addition in the lower atmosphere. 
This energy is usually assumed to come from magnetic reconnection between closed, loop-like mag- 
netic flux systems (that are in the process of emerging, fragmenting, and being otherwise jostled by 
convection) and the open flux tubes that connect to the solar wind. These models tend to explain 
the differences between fast and slow solar wind as a result of qualitatively different rates of flux 
emergenc e, reconnection, and coronal heating at the basal footpo i nts of different region s on the Sun 
(see, e.g.,'Axford and McKenzie', 'l992', 1997; 'pisk et al'., 1999; Rvutova et al!. '200l': Markovskii 
and Hollweg, 20 02. .2004: .Fisk. .2003.: .Schwadron and McComas.. 2003,;, Woo et aL.2004; Fisk and 
Zurbuchen, 200^. 

It is notable that both the WTD and RLO models have recently passed some basic "tests" of comparison 
with observations. Both kinds of model have been shown to be able to produce fast (v > 600 km/s), low- 
density wind from coronal holes and slow (v < 400 km/s), high-density wind from streamers rooted in quiet 
regions. Both kinds of model also seem able to reproduce the observed in situ trends of how frozen-in 
charge states and the FIP effect vary between fast and slow wind streams. 

The fact that both sets of ideas described above seem to mutually succeed at explaining the fast/slow 
solar wind could imply that a combination of both ideas would work best. However, it may also imply that 
the existing models do not yet contain the full range of physical processes - and that once these are included, 
one or the other may perform noticeably better than the other. It also may imply that the comparisons with 
observations have not yet been comprehensive enough to allow the true differences between the WTD and 
RLO ideas to be revealed. 
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Several recent observations have pointed to the importance of understanding the relationships and dis- 
tinctions between the WTD and RLO models. The impulsive polar jets discussed in Section [4] may be 



evidence that that magnetic reconnection drives some fraction of the fast solar wind (see also FislJ.l2005l 



Moreno-Insertis et al. 



2008tlPariat et a/.Ll2009i) . Also, direct observations of Alfven waves above the solar 



iimD maicate tne nigniy mterrmttent nature or now Kmetic energy is aistnoutea m spicules, loop s, ana tne 
open-field corona dPe Pontieu et ah . 2007 : Tomczyk et al. . 2007; Tomczyk and McIntoshL 2009 h. Spectro- 
scopic observations of blueshifts in the chromospheric network have long been interpreted as the launching 
points of solar wind streams, but it remai ns unclear how nanoflare-like events or loop-openings contribut e 



to the interpretation of these diagnostics dHe et am2007l : lAschwanden et a/.Ll2007tlMcIntosh et am2007h 
Even out in the in situ solar wind - far abo ve the roiling "furnace" of flux emergence at t he Sun - there 



remains evidence for ongoing reconnection ( Gosling et al. . 2005 : Goshng and Szabol 2008 ). There is also 
evidence that the dominant range of turbulence timescales measured in interplan etary space (i.e., tens of 
minutes to hours) is related to the timescale of flux cancellation in the low corona dHollwed . ll990Ll2006h . 

Determining whether the WTD or RLO paradigm - or some combination of the two - is the dominant 
cause of global solar wind variability is a key prerequisite to building physicall y realistic predictiv e models 
of th e heliosphere. Manv of the widelv-ap plied gl obal modeling codes (e.g.. Rilev et al. ^ 200 ll: Roussev 



eta 



20031 : Tothef a/.Ll2005l : Usmanov and GoldsteinL 12006 ; Feng et al. .L l2007h continue to utilize relatively 
simple empirical prescriptions for coronal heating in the energy conservation equation. Improving the 
identification and characterization of the key physical processes will provide a clear pathway for inserting 
more physically realistic coronal heating "modules" into three-dimensional MHD codes. 



5.2 The Wave/Turbulence-Driven (WTD) solar wind idea 

There has been substantial work over the past few decades devoted to exploring the idea that the plasma 
heating and wind acceleration along open flux tubes may be explained as a result of wave damping and tur- 
bulent cascade. No matter the relative importance of reconnections and loop-openings in the low corona, we 
do know that waves and turbulent motions are present everywhere from the photosphere to the heliosphere, 
and it is important to determine how they affect the mean state of the plasma. A review of the observational 
evidence for waves and turbulence in the solar wind is beyond the scope o f this paper, but severa l recen t 
reviews of the remote-sensing and in situ data include Tu and Marsch ('l995'),'Mullan and YakovS (Il995h . 
Goldstein et al. (1997), Roberts (2000), Bastian (2001), and Cranmer (2002a, 2004a, 2007i). Alfliough fliis 
subsection mainly describes recent work by the author, these resu l ts would not hav e been possible without 
earlier work on wave/turbulent heating by, e.g., 'Coleman' ('1968'), 'Hollweg' ('1986'), Hollweg a nd JohnsonI 
(fl9"88) . Isenberg (1990) . Li et al.. (1999) . Matthaeus et at (1999) . Dmitruk et al. (2001. ,2002.). and many 
others. 

Cranmer et al. described a set of models in which the time-steady plasma properties along a one- 

dimensional magnetic flux tube are determined. These model flux tubes are rooted in the solar photosphere 
and are extended into interplanetary space. The numerical code developed in that work, called ZEPHYR, 
solves the one-fluid equations of mass, momentum, and energy conservation simultaneously with transport 
equations for Alfvenic and acoustic wave energy. ZEPHYR is the first code capable of producing self- 
consistent solutions for the photosphere, chromosphere, corona, and solar wind that combine: (1) shock 
heating driven by an empirically guided acoustic wave spectrum, (2) extended heating from Alfven waves 
that have been partially reflected, then damped by anisotropic turbulent cascade, and (3) wind acceleration 
from gradients of gas pressure, acoustic wave pressure, and Alfven wave pressure. 

The only input "free parameters" to ZEPHYR are the photospheric lower boundary conditions for the 
waves and the radial dependence of the background magnetic field along the flux tube. The majority of 
heating in these mode l s comes from the turbulent dissipation of partially reflected Alfven w aves (see also 
Matthaeus et a/.LIl999l : lDmitruk et aZll2002l : IVerdini and 120071: [Chandran ef"aZll2009a b. Photospheric 
measurements of the horizontal motions of strong-field intergranular flux concentrations (i.e., G-band bright 
points) were used to constrain the Alfven wave power spectrum at the lower boundary. This empirically 
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determined power spectrum is dominated by wave periods of order 5-10 minutes. It is important to note, 
however, that radio and in situ measurements find that most of the fluctuation power in the solar wind 
is at lower frequencies (i.e., periods of hours). We still do not yet know (1) if the shape of the power 
spectrum evolves significantly between the lower solar atmosphere and interplanetary space, or (2) if some 
low-frequency power is missed by the existing measurements of G-band bright point motions. In any case, 
as seen below, the resulting wave reflection and turbulent dissipation that comes from just the 5-10 minute 
periods appear to be sufficient to explain the observed levels of coronal heating and solar wind acceleration. 

Non-WKB wave transport e quations were solved to determine the degree of linear reflection at heights 
above the photospheric base (see lCranmer and van Balle^ooiienll2005h . The resulting values of the Elsasser 
amplitudes Z±, which denote the energy contained in upw ard (Z-) and downwar d (Z+) propagating waves, 
were then used to constrain the energy flux in the cascade. ICranmer et al.\ (120070 used a phenomenological 
form for the damping rate that has evolved from studies of Reduced MHD and comparisons with numerical 
simulations. The resulting heating rate (in units of erg s"' cm"-') is given by 



Q = P 



1 



1 + [feddyAi-ef]"/ 4Lj_ 



ZiZ+ + ZlZ- 



(1) 



whe re p is the mass density and L± is an effective perpendicular correlation length of the turbulence (see 



e.g..lHossain et a/.lll995tlZhou and Matthaeusill990l:lBreech et a/.ll2008l:|Podesta and Bhattachari eel 120091: 



Beresnyak and Lazariar , 2009 ). ICra nmer et al! ('2007*) used a standard assumption that Lj^ scales with 
the cross-sectional width of the flux tube (HoUweg. 1986). The term in parentheses above is an efficiency 
factor that accounts for situations i n which the cascade does not h ave time to develop before the waves 
or the wind carry away the energy ("Pmitruk a nd Matthaeus , boos *). In open field regions, the cascade is 
"quenched" when the nonlinear eddy time scale fgddy becomes much longer than the macroscopic wave 
reflection time scale fief. In closed fi eld regions, the correction factor may behav e in an opposite sense as it 
does for open field regions (see, e.g., Gom ez et al. [ I2OOOI : 'RaDpazzo et a/.L l2008h . In most of the solar wind 
models, though. iCranmer eLgl-. (2007) used n = 1 in equation ([TJ based on analytic and numerical results 



'-< ■ • p-1 pj nil m ' _j |^^_^^— ± J 

dDobrowolnv et alx. Il980t lOughton et alx. 120061) . but they also tried n = 2 to explore a stronger form of the 



quenching effect. 

Figure|9]summarizes the results of varying the magnetic field properties while keeping the lower bound- 
ary conditions fixed. For a single choice for the photospheric wave properties, the models produced a real- 
istic range of slow and fast solar wind conditions. A two-dimensional model of coronal holes and streamers 
at solar minimum reproduces the latitudinal bifurcation of slow and fast streams seen by Ulysses. An active- 
region-like enhancement of the magnetic field strength in the low corona generates a high mass flux and a 
slow win d speed, in agreement with observations of open field lines connected with active regions (see also 
Wang et al .. 2009). As predicted by earlier studies, a larger coronal "expansion factor" naturally gives rise 
to a slower and denser wind, higher temperature at the coronal base, and lower-amplitude Alfven waves at 
1 AU. 

In these models, the radial gradient of the Alfven speed affect s where the wa ves are reflected and 
damped, and thus whether energy is deposited below or above thelParkei (Il958ah critical poin t. Early 
studies of solar wind energetics (e.g.. lLeer and HolzeA 1 19801 : |PneumanLll980l : lLeer et al.[ Il982h showed 
that if there is substantial heating below the critical point, its primary impact is to "puff up" the hydrostatic 
scale height, drawing more particles into the accelerating wind and thus producing a slower and more 
massive wind. If most of the heating occurs at or above the critical point, the subsonic atmosphere is 
relatively unaffected, and the local increase in energy flux has nowhere else to go but into the kinetic energy 
of the wind (leading to a faster and less dense outflow). The ZEPHYR results shown in Figure |9] display 
this kind of dichotomy because the superradial expansion creates a much higher critical point over the 
equatorial regions than over the p o les. Additional studies of how and where the mass flux an d wind speed 
are de termi ned include IWithbroe (1988), Hansteen and Leei (1995), Hansteen et al. (1997), Janse et al. 
(Hoot'), and' Wang ef a/.l(l2009h . 

Perhaps more surprisingly, varying the coronal expansion factor in the models shown in Figure |9] also 
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Figure 9: Summary of j Cranmer et al. 1 (l20Q7h models: (a) The adopted solar-minimum field geometry of 
iBanaszkiewicz et al. (1998'), with radii of wave-modified critical points marked by symbols, (b) Latitudinal 
dependence of wind speed at ^2 AU for models with n = 1 (multi-color curve) a nd » = 2 (brown curve) , 
compared with data from the first Ulysses polar pass in 1994-1995 (black curve: [Goldstein et al l fl996l) . 
(c) T{r) for polar coronal hole (red solid curve), streamer edge (blue dashed curve), and strong-field active 
region (black dotted curve) models. 



produces correlative trends that are in good agreement with in situ measurements o f commonly measured 
ion charge state ratios (e.g., O^VO^""") and FlP-sensitive abundance ratios (e.g., Fe/O). iTranmer etal. (2007) 
showed that the slowest solar wind streams - associated with active-region fields at the base - can pro- 
duce a factor of ~30 larger frozen-in ionization-state ratio of O^'^/O*''^ than high-speed streams from polar 
coronal holes, despite the fact that the temperature at 1 AU is lower in slow streams than in fast streams. 
Furthermore, when elemental fractionation is rnodeled using a theory based on preferential wave-pressure 
acceleration (lLamina . 12004 , 120091: 1 iBrvans et all . l2009h . the slow wind streams exhibit a substantial relative 
buildup of elements with low FIP with respect to the high-speed streams. Although the WTD models uti- 
lize identical photospheric lower boundary conditions for all of the flux tubes, the self-consistent solutions 
for the upper chromosphere, transition region, and low corona are qualitatively different. Feedback from 
larger heights (i.e., from variations in the flux tube expansion rate and the resulting heating rate) extends 
downward to create these differences. 

Another empirical "marker" of hehospheric stream structure is the proton specific entropy, or entropy 
per proton, which is often approximated as b eing proportional to ln(r'n/«^~^), where 7 « 1 .5 is an empirical 
adiabatic index for solar wind protons (e.g.. lBurlaga et al. iPagel et al l l2004 . When measured in 

regions of the (non-CME) heliosphere where corotating interaction regions have not yet formed shocks, this 
quantity is seen to clearly distinguish slow wind streams from fast wind streams. Figure [TO] shows how the 
specific entropy is positively correlated with wind speed, both in measurements made by the Solar Wind 
Elec tron Proton Alpha Monitor (SWEPAM) instrument on ACE dMcComas et al. ^ 1998h and in the Cranmer 



et al. (l2007h ZEPHYR models discussed above. Each model data point was computed independently of the 
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Figure 10: Solar wind specific entropy plotted as a function of solar wind speed, computed for both the 
ZEPHYR models at 1 AU (black symbols, curve) and from ACE/SWEPAM data (blue points). 



others. The models had identical lower boundary conditions at the photosphere, and they differed from 
one another only by having a different radial dependence of the magnetic field. Because entropy should 
be conserved in the absence of significant small-scale dissipation, the quantity that is measured at 1 AU 
may be a long-distance proxy for the near-Sun locations of strong coronal heating. In other words, the 
comparison of measured and modeled solar wind entropy variations may be a key way to discriminate 
between competing explanations of solar wind acceleration. 

Although equation ^ describes the plasma heating rate in terms of the local properties of MHD turbu- 
lence, it is also possible to see that this expression gives a heating rate proportional to the mean magnetic 
flux density at the coronal base. As illustrated above in FigureS] the mean field strength in the low corona 
is determined by both the photospheric field strength in the intergranular bright points and the total number 
of bright points that eventually merge their fields together in the low corona. The field strength at this merg- 
ing height can thus be estimated as B w where sa 1500 G is the (nearly universal) photospheric 
bright-point field strength and /* is the area filling factor of bright points in the photosphere. The latter 
quantity appears to vary by more than an order of magnitude in different regions on the Sun, from about 
0.002 (at low latitudes at solar minimum) to '^O.l (in active regions). If the region s belo w the merging 
height can be tr eated using approximations from "thin flux tube theorv" (e.g..lSpruiti 1981,: Cranmer and 



van Ballegooiien. 120051) . then it is possible to express each term in equation ([T]l as a function of /* and the 
photospheric properties. For example, B oc p'/^ applies to thin flux tubes in pressure equilibrium, and thus 
p at the merging height can be estimated as /^p* (where is the photospheric density). For Alfven waves 
at low heights, Z± oc p"'/"*, and so Z± at the merging height scales like Z±^,/ fl^^ . Also, we assumed that 
L±_ oc B"^/^. If the quenching factor in parentheses in equation ([T]i is neglected, then 

Equivalently, equation (|2]l implies that ~ B/B^,, and thus that the heating in the low corona scales 

directly with the mean magnetic field strength there. In a more highly structured field, the latter is equivalent 
to the magnetic "flux density" averaged over a given region. Observational evidence for such a linear scaling 
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has been found for both a variety of s olar regions and other stars as well (s ee, e.g.. lPevtsov et a/.L 12003 
Schwadron et a/.L 120061: [Suzukill2006l : lKoiima et aZl 120071: IPinto et a/.Ll2009l) . 



5.3 The Reconnection/Loop-Opening (RLO) solar wind idea 

It is clear from observations of the Sun's highlv dynamical "magnetic carpet' ' dSchriiver ef a/.lll997l: Ti- 



tle and Schrijver, .1998; Hag enaar et ai, 1999 ) that much of coronal heating is driven by the continuous 
interplay between the emergence, separation, merging, and cancellation of small-scale magnetic elements. 
Reco nnection seems to be the most likely channel for the injected magnetic energy to be converted to heat 
(e.g.. iPriest and Forbesi, ,2000) . Only a small fraction of the pho tospheric magnetic flux is in the form of 
open flux tubes connected to the heliosphere dClose et a/.ll2003h . Thus, the idea has arisen that the dom- 
inant source of energ y for open flux tubes is a series of stochastic reconnection events between the open 
and closed fields (e.g Jpisk et a/.l [l999: Rvutova et a/.lEo Ol: Fisk, 2003: Schwadron and McComas. 20031: 



Feldman et aZl 120051: ISchwadron et g/1 120061: ISchwadron and Mc Comasll2008l:lFisk and Zhaoll2009h . 



The natural appeal of the RL O idea is evident fr om the fact that open flux tubes are always rooted in 



the vicinity of closed loops (e.g.. lDowdv et a/.L 11986) and that all layers of the solar atmosphere seem to 



be in continual motion with a wide range of timescales. In fact, observed correlations between the lengths 
of cl osed loops in various region s, the electro n temperature in the low corona, and the wind speed at 1 AU 
dFeldman et a/.. 1999: Gloeckler et al. . 2003) are highly suggestive of a net transfer of Poynting flux from 
the loops to the open-field regions that may be key to understanding the macroscopic structure of the solar 
wind. The proposed RLO reconnection events may also be useful in generating energetic particles and 
cross-field diffusive transport throughout the heliosphere (e.g.. lFisk and Sch wadron. 2001). 

Testing the RLO idea using theoretical models seems to be more difficult than testing the WTD idea 
because of the complex multi-scale nature of magnetic reconnection. It can be argued that one needs 
to create fully three-dimensional models of the coronal magnetic field (arising from multiple magnetic 
elements on the surface) to truly assess the full range of closed/open flux interactions. The idea of modeling 
the coronal field via a collection of discrete magnetic sources (referred to in various contexts as "magneto- 
chemistry," "tectonics," or "ma gnetic charge topology") has been used extensively to study the evolution 
of the closed-field corona (e.g., Longcor)eL 19961: Schriiver et al. . 19971: Lon^cope and Kankelborg . 119991: 



Sturrock et a/ .'. 1999: Priest et a/.. 2002: Beveridge et a/., '2003': 'Barnes et a/., 2005: Parnell, 2007: Ng and 



Bhattachariee, 2008) but applications to open fields and the solar wind have been rarer (see, however, Fisd 



2005:Tuetal. 



20051) . 



In order to develop the RLO paradigm to the point where it can be tested more quantitatively, several 
key questions remain to be answered. For example, how much magnetic flux actually opens up in the 
magnetic carpet? Also, what is the time and space distribution of reconnection-driven energy addition into 
the (transiently) open flux tubes? Lastly, how is the reconnection energy distributed into various forms (e.g., 
bulk kinetic energy in "jets," thermal energy, waves, turbulence, and energetic particles) that each affect the 
accelerating solar wind in different ways? Combinations of simulations, analytic scaling relations, and 
observations are needed to make further progress. 



5.4 Kinetic microphysics 

The theoretical models discussed in the previous two subsections mainly involved a "one-fluid" or MHD 
approach to the coronal heating and solar wind acceleration. However, at large heights in coronal holes, 
the collisionless divergence of plasma parameters for protons, electrons, and heavy ions allows the multi- 
fluid kinetic processes to be distinguished in a more definitive way. The UVCS measurements of strong 
O'^'' preferential heating, preferential acceleration, and temperature anisotropy have spurred a great deal of 
theoretical work in this direction (see reviews by Holl weg and Isenberg . 2002; Cranmer . 2002a; Marsch, 
[2005. .2006; .Kohl et a/.>,2006i) . Specifically, the observed ordering of Ti^Tp > Tg and the existence of 
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Reconnection events in the low corona generate: 

Ion cyclotron resonant Alfven waves Tu & Marsch 1997; Cranmer 2000, 2001 

Electron beams ion cyclotron waves Markovskii & Hollweg 2002, 2004 

Fast collisionless shocks Lee & Wu 2000 



MHD turbulence in the extended corona generates: 

Ion cyclotron waves ("parallel cascade?"), with: 

^ Alfven and fast-mode nonlinear coupling Chandran 2005 

^ Three- wave (ion-sound/parametric) coupling Yoon & Fang 2008 

^ Fermi-like diffusion between inward/outward wave resonances Isenberg 2001 

Kinetic Alfven waves ("perpendicular cascade"), with: 

^ Shear instabilities ion cyclotron waves Markovskii et al. 2006 

^ Nonlinear wave-particle resonances Voitenko & Goossens 2003, 2004 

^ Debye-scale electron holes Matthaeus et al. 2003; Cranmer & van Balle. 2003 

Oblique MHD waves (high fc||, high k±) Li & Habbal 2001 

Current sheets coherent Fermi acceleration Dmitruk et al. 2004 

Transverse density gradients drift currents Markovskii 2001; Zhang 2003 



Low-frequency Alfven waves in the corona directly undergo: 

Polarization drift lower-hybrid waves Singh & Khazanov 2004, 2007 

Nonresonant stochastic heating Lu & Li 2007; Wu & Yoon 2007 

Stochastic heating at fractional eye. resonance Chen et al. 2001; Guo et al. 2008 



Heavy ion velocity filtration Pierrard & Lamy 2003; Pierrard et al. 2004 



Figure 1 1 : Tabular outline of suggested physical processes for preferentially heating and accelerating minor 
ions in coronal holes. 



anisotropies of the form T±_ > T\\ in coronal holes led to a resurgence of interest in models of ion cyclotron 
resonance. 

The ion cyclotron heating mechanism is a classical resonance between left-hand polarized Alfven waves 
and the Larmor gyrations of positive ions around the background magnetic field. If the wave frequency and 
the natural ion gyrofrequency are equal, then in the rest frame of the ion the oscillating electric and magnetic 
fields of the wave are no longer felt by the ion to be oscillating. The ion in such a frame senses a constant 
DC electric field, and it can secularly gain or lose energy depending on the relative phase between the ion's 
velocity vector and the electric field direction. In a wave field with random phases, an ion will undergo a 
random walk in energy. Thus, on average the ions can be considered to "diffus e" into faster ( i.e., wider) 
Larmor orbits with larger perpendicu lar energy (seelko wlands et al., 1966; Galins kv and Shevch enko. 200o| 
Esenbergil200ltlCranmer[|200itllsenberg and VasquezLl2009l) . 

In the actual solar corona, however, it is not likely that the situation is as straightforward as summarized 
above. Instead of a population of pre-existing, linear cyclotron waves that are dissipated, there may be 
a rich variety of nonlinear plasma mechanisms at play. The observed ion heating is likely to be just the 
final stage of a multi-step process of energy conversion between waves, turbulent motions, reconnection 
structures, and various kinds of distortions in the particle velocity distributions. Figure [TT] surveys the field 
of suggested possibilities, and the remainder of this section discusses these ideas in more detail. 

One potential obstacle to the idea of ion cyclotron heating is that the required gyroresonant wave fre- 
quencies in the corona are of order 10^ to 10"* Hz, whereas the dominant frequencies of Alfven waves 
believed to be emitted by th e Sun are thought to be much lower (i.e., less than 0.01 Hz, corresponding to 
periods of minutes to hours). Axford and McKenzie ( 1992 !) suggested that the right kinds of high-frequency 
waves may be generated in small-scale reconnection events in the chaotic "furnace" of the supergranular 
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network. These waves could propagate upwards in height - and downwards in magnetic field strength - 
until they reached a location where they becam e cyclotron resonant with the local ions, an d thus would 
damp rapidly to provide the ion heating (see also lSchwartz et a/.Ml981uTu and Marschll 19971) . 

The a bove scenario of "basal generation" of ion cyclotron waves has been called into question for several 
reasons. ICranmed (l2000l 120011) argued that the passive sweeping of a pre-existing fluctuation spectrum 
would involve ions with low gyrofrequencies (i.e., small ratios of charge to mass; « 0.1 to 0.2 in 

proton units) encountering waves of a given frequency at lower heights than the ions that have been observed 
to exhibit preferential heating, like O"""^ {qi/mi = 0.31) and Mg"^^ iqi/nii = 0.37). Thus, the resonances of 
many minor ion species may be strong enough to damp out a base- generated spectrum of waves before 
they can become resonant with the observed species. Furthermore, iHoUwe 

3 feooo) found that a base- 
generated spectrum of ion cyclotron waves would exhibit a very different appearance in interplanetary 
radio scintillations than the observed radio data. There remains some uncertainty about these criticisms of 
a basal spectrum of ion cyclotron waves, and defini tive conclusions cannot yet be made (see discussions in 



Tu and Marschll200ll:lHollweg and Isenbergll2002h . 



There are several other interesting consequences of the lAxford and McKenzid (119921) idea of rapid re- 
connection events at the coronal base. It is possible, for example, that such microflaring activity would give 
rise to intermittent bursts of parallel electron beams that propagate up into the extended corona. Sufficiently 
strong beams may be unstable to the growth of wave power at the ion cvclotron frequencies (Markovskii and 
Holl weg. .2002, 2004; Voitenko and Goossens, .2002) , and these waves would then go on to heat the ions. 
Also, iLee and Wul (l2000r) suggested that small-scale reconnection events could produce fast collisionless 
shocks in the extended corona. For shocks sufficiently thin and strong (i.e., with a bulk velocity jump of at 
least '^0.3 times the Alfven speed), ions that cross from one side of the shock to the other remain "nonde- 
flected" by the rapid change in direct ion of the magnetic fiel d. Thus, they can convert some of their parallel 
motion into perpendicular gyration. Mancuso et al. (I2OO2') suggested this mechanism may be applied to 
understanding UVCS measurements of ion heating in shocks associated with CMEs. However, it is unclear 
to what extent the open magn etic regions in coronal holes ar e filled with sufficiently strong shocks to enable 
this process to occur (see also Hollwe g and IsenbergL 2002). 

In contrast to the ideas of base-generation of ion cyclotron waves, there have been several proposed 
mechanisms for "gradual generation" of these waves over a range of distances in the corona and solar 
wind. A natural way to produce such an extended source of fluctuations is MHD turbulent cascade, which 
continually transports power at large scales to small scales via the stochastic shredding of transient eddies . 
A strong turbulent cascade i s certainly present in interplanetary space (see reviews by ITu and Marsch . 



I995I: iGoldstein et al\. \\9rh . It is well known, though, that in both numerical simulations and analytic 



descriptions of Alfven-wave turbulence (with a strong background "guide field" like in the corona) the 
cascade from large to small length scales (i.e., from small to large wavenumbers) occurs most efficiently 
for modes that do not increase in frequency. In other words, the cascade acts most rapidly to inc rease the 
perpendicular wavenumber ^ 1 while leaving the parallel wavenumber h\ largely unchanged ( e.g., StraussL 



perpenaicuiar wavenumber / c 1 wniie leavmg tne parallel wavenumber /c|| largely unc 
19761: IShebaUn et al\. [l98l lOoldreich and Sridhail [l995l : ICho et aZi l200a loughti 



:on et al 



This 

type of cascade is expected to generate so-called kinetic Alfven waves (KAWs) with k±^ 3> k\\, but not ion 
cyclotron waves. 

Under typical "low plasma beta" conditions in the corona and fast solar wind , the li near dissipation of 
KAWs would lead to the preferential parallel heating of electrons (Leamon et al. '1999^, Cranmer and van 
Ballegooijen, 2003; Gm^ and Bprpvsky, 2008). This is essentially the opposite of what has been observed 
with UVCS. However, there have been several suggestions for more complex (nonlinear or multi-step) 
processes that may be responsible for ions to receive perpendicular heating from KAW-type fluctuations. 

1. Markovskii et al. (l2006h discussed how a perpendicular turbulent cascade produces increasingly 
strong shear motions transverse to the magnetic field, and that this shear may eventually be un- 
stable to the generation o f cyclotron resonant waves that can in turn heat protons and ions (see also 



Mikhailenko et al 



2008V This effect may also produce a steepening in the power spectrum of the 
magnetic field fluctuations that agrees with the observed "dissipation range" tSmith et aL,2006) . 
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Voitenko and Goossens (120031120041) suggested that high-fcj^ KAWs with sufficiently large amplitudes 
could begin to exhibit nonUnear resonance effects ("demagnetized wave phases") leading to rapid ion 
perpendicular heating. There are definite thresholds in the KAW amplitude that must be exceeded for 
these mechanisms to be initiated, and it is unclear whether the actua l coronal turbu l ent spectrum has 
enough power in the relevant regions of wavenumber space (see also Wu and YangL 20061 2007 ). 



The damping of low-frequency KAWs may give rise to substantial parallel electron acceleration. If 
the resulting electron velocity distributions were sufficiently beamed, they could become unstable to 
the generation of parallel Langmuir waves. In turn, the evolved Langmuir wave trains may exhibit a 
periodic electric potential-well structure in which some of the beam electrons can become trapped. 
Adjacent potential wells can t hen merge with on e another to form i solate d "electron phase space 
holes" of saturated p otential. lErgun et al. I (Il999h. iMatthaeus et al.l I (l2003l). and Cranmer and van 
Ballegooijen (2003) described how these tiny (Debye-scale) electrostatic structures can heat ions 
perpendicularly via Coulomb-Uke quasi-coUisions. 

Obliquely propagating MHD waves with large perpendicular and parallel wavenumbers - including 
KAWs and fast-mode waves - can interact resonantly with positive ions via channels that are not 
available when either or k± are small. Li a nd Habbal (2001) found that oblique fast-mode waves 
with large wavenumbers may be even more e fficient than Alfven waves at heating ions under coronal 
conditions. iHollweg and Markovskiil (2002) discussed how the higher-order cyclotron resonances 
become available to obliquely propagating waves with large wavenumbers, and how these can lead 
to stochastic velocity-space diffusion for ions. 

On the smallest spatial scales, the plasma in numerical simulations of MHD turbulence is seen to 
develop into a collection of narrow current sheets undergoing oblique m agnetic reconnection (i.e., 
with the strong "guide field" remaining relatively unchanged). Dmitruk et al. I (12004.) performed 
test-particle simulations in a turbulent plasma and found that protons can become perpendicularly 
accelerated around the guide field beca use of co herent forcing from the perturbed fields associated 
with the current sheets (see also iParashar et al.l 12009,) . It remains to be seen whether this process 
could lead to more than mass-proportional energization for minor ions. 

If the plasma contains sufficiently small-scale density gradients transverse to the magnetic field 
(V±p), then drift currents can be excited that are unstable to the generation of high-frequen cy waves 
jMarkovskiii I2OOII: IZhangl 120031: IVranies and Poedtsl I2OO8I: iMecheri and Marschl I2OO8I) . These 
instabiUties depend on both the amplitudes and scale lengths of \/'±p. To measure the latter, it is im- 
porta nt to take in t o account both re mote-sensing measurements of coronal density inhomogeneities 
(e.g., l2006l:IPasachoff ef ad l2007p and constraints from radio scin tillation power spectra at 
larger distances dBastianl 2001 ; Spangleii 2002 ; Harmon and Colesl 2005 ). 



Despite the fact that theory predicts a predominantly perpendicular cascade, there is some evidence 
that the turbulent fluctu ations in the solar wind have some energy that extends up to large kn values in a 



power-law tail (see, e.g., Bieber et al. , 19961 : Dasso et al. , 2005 ; MacBride et al. , 2008 ). Whether or not this 



means that true "parallel cascade" occurs in the corona and solar wind is still not known. However, some 
progress has been made using a phenomenological approach to modeUng the cascade as a combinatio n of 
advection and diffusion in wavenumber space. In the model of Cranmer and van Ballegooijen ( 2003 h. the 
relative strengths of perpendicular advection and diffusion determine the slope of the power-law spectrum 
in and thus they specify the amount of wave energy that is available at the ion cyclotron frequencies (see 



in /C|i , ana tnus tney speciiy tne amount 01 wave energy tnat is avaiiaoie at 
also lCranmer et adll999aLlLandi and Cranmeill2009l : lJiang et al. ■ I2OO9I ) 



There have also been proposals for additional mechanisms that could allow a parallel cascade to occur 
in the corona and solar wind. No nlinear couplings between the dominant A lfven waves and other modes 
such as fast magnetosonic waves ('Chandran'.'2005':^ Luo and M elrosel '2006*) and ion-acoustic waves (Yoon 
and Fang. , 2008.) have the potential to enhance the wave power at high frequencies. It has also been known 
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for some time that nonlinear coupling between Alfven and fast-mode waves may help explain why the 
measured in situ magnetic f i eld magnitude |Bj remains roughly constant while its direction varies strongly 



measured in situ magnetic neid magnituae iBj rernains rougniy c 
dBarnes and Hollwe 3. ll974l:IVasauez and Hollw"e3.ll996lll998h . 

If MHD waves have sufficiently large amplitudes, they may undergo nonlinear wave steepening, which 
leads to density variations as well as oscillations in the parallel components of the velo city and magnetic 
field. These may ge nerate progressively smaller scales along the magnetic field (e.g., iMedvedevL I2OOOI: 



Suzuki et all 120071). There is also a "bootstrap" kind of effect for ion cyclotron wave generation that 



was discussed bv llsenberg et al.\ (1200 Ih . If some outward-propagating cyclotron waves exist, the resonant 
diffusion may act to produce proton velocity distributions that are unstable to the generation of inward- 
propagating cyclotron waves. In response, the proton distributions would become further deformed and 
thus could become unstable to the growth of both inward and outward waves. It is not yet known if this 
process could reach the point of becoming self-sustaining, but if so, it may also serve as an extended 
generation mechanism for high-^|| waves. 

In addition to the above ideas that involve large wavenumbers and kinetic effects, there have been other 
suggested physical processes that do not require high-A; resonances to be initially present in order to heat 
the ions. 

1. Particles in large-amplitude Alfven waves exhibit both E x B drift motions (i.e., their standard ve- 
locity amplitude) and a polarization drift velocity Vpoi that is smaller than the former by the ratio 
uj I il,, where w is the wave frequency and fi,- is the ion cyclotron frequency. A sufficiently large Vpoi 
can lead to cross-field currents unstable to the gener ation of high-frequency waves, and to eventual 
equipartition between Vpoi a nd the ion thermal speed ( Singh and Khazanov, 2004; Sin gh et ai , IOOtI 



Khazanov and Singhll2007r) . It is not yet known whether the effective Vpoi for the coronal fluctuation 



spectrum is large enough to provide a significant fraction of the ion thermal speeds. 

Recently there have been suggested some completely nonresonant mechanisms that de pend on the 
stochasticity of MHD turbulence to produce an effec tive increase in random ion motions (ILu and LiL 



2007l : IWu and YoonLl2007l : lBourouaine et a/.ll2008 ). Questions still remain, though, concerning the 
spatial scales over which one should refer to particle motions as "heating" versus "wave sloshing." 
This energization mechanism may be just a more chaotic form of the standard velocity amplitude tha t 
an ion feels when in the presence of a spectrum of Alfven waves (see, e.g.. lWang and Wul l2009h . 
In this case, the maximum amount of heating from this process would provide mass-proportional 
heating for minor ions and protons (i.e., Tj/Tp = mi/nip), and it is clear that the UVCS measurements 
for O"^^ show heating in excess of this amount (see Section[ 



Both numerical and analytic studies of Alfven waves show that, at sufficiently large amplitudes, there 
can be gyroresonance-like ion energizati on for sets of frequencies at specific fractions of the local 



ion cyclotron resonance frequency (e.g.. lChen et g/.lBOOlHGuo et a/.L 120081) . Like several other 



proce sses listed above, this ef fect becomes active only above certain thresholds of wave amplitude. 
Also, iMarkovskii et al. I (l2009h showed that mildly nonlinear Alfven waves - with frequencies slightly 
below the local proton gyrofrequency and power in both the upward and downward directions along 
the field - can also undergo additional modes of dissipation and proton heating that are not anticipated 
in linear ion cyclotron resonance theories. 

Finally, there has been some development of the so-called velocity filtration theory, which requires nei- 
ther direct heating nor wave damping in order to energize coronal ions. Spacecraft measurements of plasma 
velocity distributions, both in the solar wind and in planetary magnetospheres and magnetosheaths, have 
re vealed tha t "supratherma l" power-law tails are quite common. These observations led to the suggestion 
bv lScuddel^(^1992 a'■'l994^ of an an alternative to theories that demand explicit energy deposition in the 
low corona (see also Parker, 1958b; Levine, 1974). A velocity distribution having a suprathermal tail will 
become increasingly dominated by its high-energy particles at larger distances from the solar gravity well. 



25 



Thus an effective "heating" occurs as a result of particle-by-particle conservation of energy. The major 
unresolved issue is whether suprathermal tails of the required strength can be produced and maintained in 
the upper chromosphere and transition region - where Coulomb collisions are traditionally believed to be 
strong enough to rapidly drive velocity distributions toward Maxwellians. Whether the solar atmosphere 
actually plays host to strong nonthermal tails is still under deba t e, with some evidence existing in favor 
of their presence (e.g., Esser and Edgar, 2000; Ralchenko et al. . 2007 ) and other evidence against them 

996: Feldman et a/.. .2007.) . 

b|) ideas about supratherm al velocity filtration were a ppli ed only to the 



(lAnderson et a/.LI1996tlKo et al. 



The original dScuddeii Il992al 



primar y (proton and electron) coronal plasma. More recentlv. lPierrard and Lamv ( 200 3") and Pi errard et al. 



(120041) have shown that this mechanism can produce extremely high temperatures for heavy ions in the 
corona - providing they had suprathermal tails in the chromosphere. The primary quantities presented 
in these papers, however, were integrated isotropic temperatures T . No information was given about the 
predicted sense of the temperature anisotropy for the minor ions. For a collisionless exospheric model, 
there is a suspicion that a combination of several effects (e.g., the initial velocity filtration and subsequent 
magnetic moment conservation) would result in velocity distributions with T\\ ^ T^, which is not what is 
observed. 

Since it is obvious that not all of the proposed mechanisms described above (and shown in Figure fTTTi 

can be the dominant cause of the collisionless ion energization in coronal holes, there is a great need to 

"cut through the jungle" and assess the validity of each of these processes. For many of these suggested 

ideas, further theoretical development is required so that specific observational predictions can be made. 

However, there are also several types of measurement that have not been widely recognized or utilized 

as constraints on theoretical models. A prime example is the use of radio sounding (i.e., interplanetary 

scintillations and Faraday rotation) to measure the fine structure of the corona and solar wind in density, 

1 ' 1 ' 1 1 — — I 

velocity, and magnetic field strength (see, however, Hollweg and Isenberg, 2002; Spangler, 2002; Harmon 

and Coles, 2005 ; Chandran et al-, 2009b). Another example is the use of high-resolution UV sp ectral line 

profil es to probe departures from Maxwellian or bi-Maxwellian ion velocity distributions (e.g.. lCranmer . 



2001). 



6 Summary and Conclusions 

The last decade has seen significant progress toward identifying and characterizing the processes that pro- 
duce coronal holes. As remote-sensing plasma measurements have become possible in the extended solar 
corona (i.e., the region of primary acceleration of the solar wind), the traditional gap between solar physics 
and in situ space physics has become narrower However, there are still many unanswered questions: How 
and where in the solar atmosphere are the relevant waves and turbulent motions generated? Which kinds of 
fluctuation modes (i.e., linear or nonlinear; Alfven, fast, or slow; high k\\ or high k±) are most important? 
What frequencies dominate the radially evolving power spectrum? What fraction of the interplanetary solar 
wind comes from filamentary structures such as polar plumes and polar jets? Are there relatively simple 
"scaling laws" that will allow us to use only the measured properties at the solar surface to predict the 
resulting amount of coronal heating and solar wind acceleration? 

Answering the above questions involves moving forward in both the theoretical and observational di- 
rections. Section|5]described the large number of suggested physical processes for energizing the plasma in 
coronal holes. The validity of many of these processes still needs to be assessed, and their relative contri- 
butions to the heating and acceleration of the actual solar corona need to be quantified. If, at the end of this 
process, there are still a number of mutually inconsistent theories that are still viable, the only way forward 
is to determine what future measurements would best put the remaining models to the test. These activ- 
ities are ongoing with the planning of inner heliosph eric missions such as Solar Probe (McComas et a/I 



20071) and Solar Orbiter dMarsden and Fleckl 120071) . as well as next-generation ultrav iolet coronagraph 



spectroscopy missions that would follow up on the successes of mCSISOHO (see, e.g.. lCranmeiil2002bl: 
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Pmtostar CTTS WTTS ZAMS Present-day 

{-iO^yr) (-10^ yr) (-lo" yr) (-10* yr) Sun (-10^ yr) 



Figure 12: Illustration of the evolving circumstellar environment of a solar-mass star (see text), showing 
various kinds of open-field structures that may be analogous to present-day coronal holes. 



Gardner et a/.ll2003l:lKohl et a/.Ll2006h . 



The development of more physically sound models of the solar wind feeds back in many ways to a 
wider understanding of stellar outflows and star/planet evolution. Figure [12] shows some of the the early 
stages of evolution for a represe ntative solar-type star. At all ages, cool stars are inferred to exhibit sorn e 
kind of wind or jet-like outflow ( Lamers and Cassinelli , 1 1 9991: 1 Woodl |2004 lOiidell l2007l: ICranmer . 2008 
Young stars first become visible as dust-obscured cloud cores and protostars (e. g.. iLadal 1 1 985» : ,Hartman 
[2000), and these objects are often associated with bipolar, collimated jets. These outflows indicate some 
kind of transfer of energy from the accretion disk's orbital motion to torqued magnetic fields (rooted on 
the stellar surface) that relieve the buildup of angular momentum and eject plasma out the poles (e.g.. 



lie 

I 



Blandford and Payn e. 1982; van Ballegooijen, 1994). As the accretion rates decrease over time, protostars 
become visible as classical T Tauri stars (CTTS), an d there remains ample evidence for p olar outflows in 



)r polar o 

the for m of both "disk winds" and true stellar winds ('Hartig an et ali 119951 : iFerreira et al. 1 120061: ICranmeA 
2008bh . The primordial accretion disk is dissipated gradually as the star enters the weak-lined T Tauri 



star (WTTS) phase, with a likely transition to a protoplanetary dust/debris disk. Strong stellar m agnetic 
activity remains evident during these stages from, e.g.. X-rays (^ Feigelson and Montmerle . 1999h . Many 
"post T Tauri" stars, once they reach the zero-age main sequence (ZAMS), remain rapidly rotating, and 
for young ZAMS stars such as AB Dor there is evidence for a range of X-ray emitting plasma from dark 
polar spots (which probably do not correspond to open magneti c field regions like coronal holes) to huge 
"slingshot promi nences" extending over several stellar radii (e.g.. Giidel et al. ^ 2001 ^ 2003: Jardine and van 
Ballegooijen. i2005h . 

Learning about the fundamental physics responsible for solar coronal holes has relevance that reaches 
into other areas of study besides astrophysics, including plasma physics, space physics, and astronautical 
engineering. The practical benefits of improvi ng long-term predictions for the cond itions of the Earth's 
local space environment are manifold (see, e.g.. lFeynman and GabrielL 2000l : Eastwood i2008.) . In addition, 
parallel research into the expansion of the polar wind from the Earth's ionosphere has led to an improved 
understa nding of kinetic processes in plasmas on the bou ndary between collisional and collisionless con- 



ditions (.Lemaire and PierrardluOOluBarakat and Schunld,l2006i) . Finally, a growing realization that strict 



topical compartmentalization is often a hindrance to making pro gress has given rise to greater inte rest in 



interdisciplinary studies of "universal processes in heliophysics" ( Crooker , 2004 ; Davila et al. , 2009h . 
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